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INTEGRAL FIELD SPECTROSCOPY OF THE EXTENDED EMISSION-LINE REGION OF 4C 37.43* 
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ABSTRACT 

We present Gemini integral field spectroscopy and Keck II longslit spectroscopy of the extended emission- 
line region (EELR) around the quasar 4C 37.43. The velocity structure of the ionized gas is complex and cannot 
be explained globally by a simple dynamical model. The spectra from the clouds are inconsistent with shock or 
"shock + precursor" ionization models, but they are consistent with photoionization by the quasar nucleus. The 
best-fit photoionization model requires a low-metallicity (12+log(0/H) < 8.7) two-phase medium, consisting 
of a matter-bounded diffuse component with a unity filling-factor (N ~ 1 cm~ 3 , T ~ 15000 K), in which 
are embedded small, dense clouds (N ~ 400 cm" 3 , T ~ 10 4 K). The high-density clouds are transient and 
can be re-generated through compressing the diffuse medium by low-speed shocks (Vs < 100 km s" 1 ). Our 
photoionization model gives a total mass for the ionized gas of about 3 x 10 10 M Q , and the total kinetic energy 
implied by this mass and the observed velocity field is ~ 2 x 10 58 ergs. The fact that luminous EELRs are 
confined to steep-spectrum radio-loud QSOs, yet show no morphological correspondence to the radio jets, 
suggests that the driving force producing the 4C 37.43 EELR was a roughly spherical blast wave initiated by 
the production of the jet. That such a mechanism seems capable of ejecting a mass comparable to that of the 
total interstellar medium of the Milky Way suggests that "quasar-mode" feedback may indeed be an efficient 
means of regulating star formation in the early universe. 

Subject headings: quasars: individual(4C 37.43) — quasars: emission lines — galaxies: evolution — galaxies: 
ISM — galaxies: abundances 



1. INTRODUCTION 

A number of low-redshift QSOs show luminous extended 
emission-line regions (EELRs) that have characteristic di- 
mensions of a few tens of kpc, considerably larger than the 
typ ical ~ 1 kpc of the cl assical narrow-line regions (NLRs; 
see Stockton et al. 2006b for a recent brief review of EELRs). 
In a few cases, such extended emission may be simply due to 
QSO photoionization of the in situ interstellar medium (ISM) 
of the host galaxy or of nearby gas-rich dwarf galaxies, but, 
for the most luminous examples, such explanations are in- 
adequate. These luminous EELRs typically show complex 
filamentary structures that bear no close morphological rela- 
tionships either with the host galaxies or with extended radio 
structures. 

In spite of this general lack of attention on the part of 
EELRs to the structural parameters of their host galaxies, they 
quite clearly are not oblivious to t he pro pertie s of the QSOs 
themse lves. As Boroso n & Okd (119841) and iBoroson et al.l 
d 1985b found, and IStockton & MacKentvl (1 1987b confirmed. 

*BASED ON OBSERVATIONS OBTAINED AT THE GEMINI OBSER- 
VATORY, WHICH IS OPERATED BY THE ASSOCIATION OF UNIVER- 
SITIES FOR RESEARCH IN ASTRONOMY, INC., UNDER A COOPER- 
ATIVE AGREEMENT WITH THE NSF ON BEHALF OF THE GEMINI 
PARTNERSHIP: THE NATIONAL SCIENCE FOUNDATION (UNITED 
STATES), THE PARTICLE PHYSICS AND ASTRONOMY RESEARCH 
COUNCIL (UNITED KINGDOM), THE NATIONAL RESEARCH COUN- 
CIL (CANADA), CONICYT (CHILE), THE AUSTRALIAN RESEARCH 
COUNCIL (AUSTRALIA), CNPQ (BRAZIL) AND CONICET (AR- 
GENTINA). SOME OF THE DATA PRESENTED HEREIN WERE OB- 
TAINED AT THE W.M. KECK OBSERVATORY, WHICH IS OPERATED 
AS A SCIENTIFIC PARTNERSHIP AMONG THE CALIFORNIA IN- 
STITUTE OF TECHNOLOGY, THE UNIVERSITY OF CALIFORNIA 
AND THE NATIONAL AERONAUTICS AND SPACE ADMINISTRA- 
TION. THE OBSERVATORY WAS MADE POSSIBLE BY THE FINAN- 
CIAL SUPPORT OF THE W.M. KECK FOUNDATION. 
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luminous EELRs are associated virt ually exclusively with 
steep- spectrum radio-loud QSOs. IStockton & MacKentvl 
dl987b also noted a strong correlation between the strength 
of the nuclear narrow-line emission and that of the extended 
emission. This latter correlation may not be particularly sur- 
prising. However, the restriction of luminous EELRs to steep- 
spectrum quasars, together with the general lack of correspon- 
dence between the distribution of the ionized gas and the ra- 
dio structure, presents an intriguing puzzle. To add to the 
confusion, not all low-redshift, steep-spectrum quasars (even 
among those that have luminous classical NLRs) show lumi- 
nous extended emission: the fraction that do is ~l/3 to 1/2 
( IStockton & MacKentvlll987l) . 

Early guesses regar ding the origin of the extende d gas cen- 
tered on ti dal debris (IStockton & M acKentvl 1 1987b or cool- 



ing flo ws dFabian et al.lll987b . However, as ICrawford et all 
(1988) correctly pointed out, ionized gas from a tidal en- 
counter would dissipate on a time scale of < 10 6 years unless 
it were confined gravitationally or by the pressure of a sur- 
rounding medium. On the other hand, Stoc kton et al.l (12002b 
showed, from a detailed photoionization model of the EELR 
around the z = 0.37 quasar 4C 37.43, that the pressure of any 
surrounding hot gas was too low for the gas to cool in less 
that a Hubble time. This conclusion was supported by the 
lack of evidence for a general distribution of hot gas around 
3 quas ars in deep Chandra X-ray imaging (Stockt on et al.l 
l2006ab . Gravitational confinement of the amount of ionized 
gas seen in many EELRs would imply unreasonably large 
masses of neutral or otherwise invisible material in the out- 
skirts of the host galaxy. 

Given the short lifetime of unconfined ionized filaments and 
the unlikelihood of gravitational confinement, along with the 
other constraints we have mentioned, the most likely scenario 
for producing an EELR is that of a superwind (IStockton et al.l 
l2002t iFu & Stockton! 120061) . In principle, such a superwind 
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could result either from a starburst or from feedback fro m the 
quasar (e.g., iDi Matteo etaT1l2005t iHopkins et aT1l2006l). with 
some evidence favoring the latter alternative dFu & Stockton! 
2006). Such an origin would be of considerable interest, since 
it would mean that we have local examples of a process sim- 
ilar (though almost certainly not identical) to the "quasar- 
mode" feedback mechanism that may be important in initially 
establishing the observed correlation between bulge mass and 
black hole mass during galaxy formation in the early universe. 
In order to explore this possibility in more detail, we have car- 
ried out an extensive re-examination of the EELR of 4C 37.43. 

4C 37.4 3 has the most lumino us EELR among QSOs at 
z < 0.45 (IStockton & MacKentylll987l) . The continuum im- 
age of the host galaxy shows distortions and a likely tidal tail 
that indicate that a major merger is in progress, which has 
probably triggered the current episode of quasar activity. Like 
most others, this EELR is morphologically completely uncor- 
rected with both the stellar distribution in the host galax y 
and with the FR II double radio source (iMiller et alJll993l) . 
It has two main condensations of ionized gas at roughly the 
same distance (~3.5", i.e., about 18 kpc) to the east and 
west of the quasar (Fig. [T). The optical spectrum of the 
brightest region, the main east condensation (El; following 
IStockton et alj 120021). has been discussed in several papers 
before (IStocktonll 19761 11 oroson & Okd[l984T: IStockton et all 
2002). Here we present our new GMOS integral field spec- 
troscopy and DEIMOS longslit spectroscopy covering a large 
fraction of the whole EELR. Throughout we assume a flat cos- 
mological model with Hq = 70 km s _1 Mpc~\ f2 OT = 0.3, 
and fi A = 0.7. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. GMOS Integral Field Spectroscopy 

4C 37.43 was observed w ith the Integral Field Unit (IFU; 
lAtlington-Smifh et al . 2002) of the Gemini Multiobject Spec- 
trograph (GMOS: lHook et al.l2004h on the Gemini North tele- 
scope. We observed a 13" x 9" region centered on the QSO 
in the early half-night of 2006 May 23 (UT). Since the main 
purpose of these observations was to determine the global ve- 
locity field of the EELR through the [O III] A5007 line, the 
IFU was used in the full-field (two-slit) mode and dithered 
on a rectangular grid of 6" x 4". The QSO was successively 
placed at each corner of the 7" x 5" IFU field and ^O'.'S away 
from the edges, so the fields covered by the pointings over- 
lap with each other by ~1" (see Fig. [TJ. We obtained two 
720-s frames and one 156-s frame on the SE pointing 2 and 
three 720-s exposures on the other three pointings. With the 
R831/G5302 grating and a central wavelength of 658.5 nm, 
we obtained a dispersion of 0.34A per pixel, a spectral resolu- 
tion of ~ 1.3 A (58 km s _1 ) FWHM, and a wavelength range 
of 6330-6920 A. The r and RG610 filters were used to avoid 
spectral overlaps. The spectrophotometric standard star Feige 
34 was observed for flux calibration. The seeing was ~ 0'.'4 
throughout the half-night. Before running the data through the 
reduction pipeline (see the 3rd paragraph in this section), the 
exposures for each pointing position were combined using the 
IRAF task IMCOMBINE. We weighed the frames according 
to their integration times. Pixels were rejected if their values 
were 1-a off the median 3 . For all four positions, there are no 
apparent cosmic rays in the final image. The shortened expo- 

2 The third exposure was cut short due to a dome problem. 

3 a estimated from known CCD parameters. 



sure of the SE pointing did make the data a little bit shallower 
in this region than others. We corrected for this difference 
by using an exposure map while making the final mosaicked 
datacube. 

We also obtained deep integral field spectroscopy ona~ 
3"5 x 5" region about 2" S W of the QSO. A total of five 2400- 
s exposures were taken with the half-field (one-slit) mode in 
the first half-night of 2006 May 24 (UT). This configuration 
provides a wider spectral coverage but a smaller field than 
does the two-slit mode. With the B600/G5303 grating and 
a central wavelength of 641.2 nm, the dispersion and spec- 
tral resolution were approximately 0.46A per pixel and 1.8 A 
(FWHM). The wavelength range was 4250-7090 A, which 
includes emission lines from [Ne V] A3426 to [Olll] A5007. 
BD +28° 4211 was observed for flux calibration. The seeing 
was ~(y/6 during the half-night. 

The data were reduced using the Gemini IRAF package 
(Version 1.8). The data reduction pipeline (GFREDUCE) 
consists of the following standard steps: bias subtraction, 
cosmic ray rejection, spectral extraction, flat-fielding, wave- 
length calibration, sky subtraction, and flux calibration. Spec- 
tra from different exposures were assembled and resampled 
to construct individual datacubes (x,y,A) with a pixel size of 
0705 (GFCUBE). For each datacube, differential atmosphere 
refraction was corrected by shifting the image slices at each 
wavelength to keep the centroid of the quasar constant. The 
four datacubes resulting from the four IFU two-slit pointings 
were then merged to form a single datacube (the IFU2 dat- 
acube), and the five datacubes from the IFU one-slit mode 
were combined to form the IFUR datacube. Finally, these 
two datacubes were binned to 0'.'2 pixels, which is the origi- 
nal spatial sampling of the IFU fiber-lenslet system. 

Since our study focuses on the emission-line gas, it is 
desirable to remove the light of the quasar from the dat- 
acubes. Since the emission-line clouds show essentially pure 
emission-line spectra, we used the continuum on either side of 
an emission line to precisely define the PSF of the quasar. We 
then ran the two-comp onent deconvoluti on task PLUC Y (con- 
tributed task in IRAF; iHook et alii 19941) on each image slice 
across the emission line to determine the flux of the quasar 
at each wavelength. Finally the quasar component was re- 
moved by subtracting a scaled PSF (according to the fluxes 
determined by PLUCY) from each image slice. Note that the 
deconvolved images produced by PLUCY were not used to 
form the QSO-removed datacubes. 

2.2. DEIMOS Longslit Spectroscopy 

We obtained additional longslit spectroscopy using 
the DEep Imag ing Multi-Object Spectrograph (DEIMOS; 
iFaberet al.ll2003l) of the Keck II telescope on 2006 August 
28 (UT). We initially centered the 20" long slit on the QSO 
at a position angle of 19° (N to E), then offset it by 0773 so 
that it goes through two EELR clouds (Nl & SI; Fig. [D that 
are presumably assoc iated with X-ray emission (see Fig. 5 in 
Stock ton" et al.ll2006ah . The alignment accuracy of DEIMOS 
is typically ~0"1 4 . The total integration time was 3600 s. 
With a 600 groove mm -1 grating, a 0'.'9 wide slit and a cen- 
tral wavelength at 7020A, we obtained a dispersion of 0.63 
A per pixel, a resolution of 3.5 A, and a spectral range of 
4570-9590 A, covering emission lines from [Ne v] A3346 to 
[S II] A6731. The spectrophotometric standard star Wolf 1346 

4 http://www2.keck.hawaii.edu/inst/deimos/specs.html 
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FIG. 1 .— Hubble Space Telescope (HST) WFPC2 [O III] image of 4C 37.43 I Stockton et al. 2002). Overlaid are the field-of-view (FOV) of the mosaicked 
GMOS IFU2 datacube (outer solid box), the field of the GMOS IFUR datacube (inner solid box), the position of the DEIMOS slit (long dashed rectangle) and 
that of the LRIS slit (long dash dotted lines). The tick marks on the inner edge of the total IFU2 FOV show the overlapping regions of the four IFU2 pointings; 
from each corner, the rectangle defined by the farthest tick mark from the corner gives the FOV for a single pointing. The radio jet direction is described by the 
solid lines. The cross near the center marks the position of the QSO. 



was observed with a 1'.'5 wide slit at parallactic angle. The 
GG400 filter was used for all observations. The seeing was 
~ 0'/7 during the observations. We reduced the data using 
the spec2d data reduction pipeline 5 . The standard star spec- 
trum suffers second-order contamination in the spectral range 
> 8000 A; for the spectra of the EELR this contamination 
is negligible since there is essentially no continuum. How- 
ever, in order to flux calibrate the emission-line spectra, one 
needs to correct for the second-order contamination of the 
standard star spectrum. We thus obtained a DEIMOS spec- 
trum of Feige 34 taken on November 29 2006 (UT). This stan- 
dard star was observed using the same settings as ours except 
that a different filter (GG495) and a slightly different central 
wavelength were used. Since both GG400 and GG495 filters 
reach an identical transmission value redward of 7000 A, no 
correction for the transmission curves needs to be applied. We 
first derived the calibration files with STANDARD for both 
spectra, then multiplied the Feige 34 data by a constant so 
that the two matched between 7000 and 7600 A. Then, in 
the bluer region, we selected the data points from the GG400 
spectrum (Wolf 1346), while in the redder region, we selected 
the GG495 data (Feige 34). We then fitted a smooth curve to 
the selected data to form the final sensitivity function (SENS- 
FUNC), which was used in the flux calibration. 

Proper sky subtraction is often a problem beyond 7000 A 
because of strong airglow lines and fringing in the CCD. 
The spec2d pipeline attempts to defringe the spectra by di- 
viding the science spectra with a normalized flat field. This 

5 The specld pipeline was developed at UC Berkeley with support from 
NSF grant AST-0071048. 



feature significantly improves the sky subtraction result, al- 
though it is still not perfect when dealing with the strongest 
sky lines. Fortunately, the redshift of 4C 37.43 nicely places 
the important red lines (i.e., [Nil] AA6548,6584, Ha and 
[S II] AA6717,6731) in a window (8900-9300 A) between two 
strong OH bands that is essentially free of strong sky lines, 
making the subtraction of sky lines almost perfect in this re- 
gion. 

Figure[T]presents an overview of the regions covered by the 
observations discussed above. We also include a deep Keck 
I LRI S (Low-Resol ution Imaging Spec trometer; lOke et all 
1 19951) spectrum from lStockton et all (120021) in our discussion. 
That spectrum had a total integration of 3600 s and was taken 
with a 300 groove mm -1 grating and a 1" wide slit. The po- 
sition of this slit is also shown in Fig.Q] 

3. RESULTS 

3.1. Kinematics 

Figure |2] presents the velocity field of the [O III] A5007 
emission derived from the IFU2 datacube, which has a spec- 
tral range that brackets the H/3-[0 III] region. The QSO was 
r emov ed from the datacube using the method described in 
§ 12.11 We discarded spectra that show very low amplitude- 
to-noise ratios (A/N < 3), and spatially binned the remaining 
spectra with A/N < 10 to a target A/N ~ 8 us ing the Voronoi 
binning method of ICappellari & Co pin (2003). We then used 
multiple Gaussians to fit the [O III] A5007 line profile in the 
binned datacube. Because different velocity components are 
sometimes present along the same line of sight, the velocity 
field is displayed in three velocity bands from —620 to +250 
km s _1 (relative to zq = 0.37120, the redshift of the quasar 
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FIG. 2.— Velocity field of 4C 37.43 EELR derived from the [O III] A5007 region of the QSO-subtracted IFU2 datacube. The three columns, from left to right, 
are line intensity, radial velocity (relative to that of the nuclear NLR) and velocity dispersion maps. To separate different clouds that are present along the same 
line of sight, the velocity field is shown in three separate velocity bands (—620 to —170 km s , —290 to +80 km s _1 and —70 to +250 km s — 1 from top to 
bottom), and the line intensity and velocity dispersion m aps r efer to the same clouds. Pixels are tf.'2 squares. The crosses indicate the position of the quasar 
before it was removed from the datacube (for details see § 12. U . Contours are from the line intensity maps. 



NLR as determined from its [O III] AA4959,5007 lines). The 
velocity dispersion (a) measurements were corrected for the 
do = 24.6 km s _1 instrumental resolution. Monte-Carlo 
simulations were performed to determine the uncertainties in 
the measured line parameters. For A/N = 8 and an intrinsic 
a = 50 km s _1 (130 km s -1 ), the l-cr errors on velocities 
(V) and velocity dispersions (a) are both about 4 km s _1 (6 
km s _1 ), and the errors on fluxes are about 6% (4%). 

Our velocity maps are generally consistent with those from 
previous i ntegral-field spectroscopy and image-sliced spec- 
troscopy (iDurret et all [19941 ICrawford & Vanderriestl 120001: 
IStockton et alj |2002). but they offer a much clearer view of 
the complexities in the velocity structure of this region. Over- 
all, these maps show that: (1) the clouds comprising the EELR 
are certainly not in a coherent rotat ion about the QSO, as al- 
ready pointed out bv lStocktonld 19761) with very limited veloc- 
ity information; (2) the majority of the clouds are blueshifted 
relative to the QSO, and the highest blueshifted velocity is 
about —620 km s — 1 , while the highest redshifted velocity is 
only +250 km s -1 ; (3) the velocity dispersions for the most 
part are between 50 and 130 km s _1 (or 120 to 310 km s _1 
FWHM), much higher than the sound speed in a hydrogen 
plasma at 10 4 K, c s ~ 17 km s -1 ; (4) there is no obvious ev- 
idence for jet-cloud interactions. As seen in Fig.Q] the radio 
jet just misses the two major condensations, and no signifi- 



cant increase in velocity dispersion is observed where the jet 
crosses any visible part of a cloud. 

As shown in the —620 to —170 km s _1 panels, the two 
major condensations to the east (El) and west of the QSO 
are both blueshifted to about —250 km s — 1 . El is resolved 
in our velocity maps, thanks to the Of! 2 sampling of the IFU 
and the Of! 4 seeing. The southern half of El shows a velocity 
increase from — 180 to —350 km s _1 along the W-E direction, 
and also an increase in a from 50 to 110 km s _1 along the 
same direction. This kinetic structu re suggests that El is not a 
coherent body (first suspected by Durret et al. 1994); instead, 
it is composed of sub-clouds, as clearly seen in the HST [O III] 
image (see the low-contrast insets to Fig. 3 of IStockton et al.l 
l2002h . 

Overall, the velocity field of the 4 C 37.43 EELR is sim ilar 
to that of the EELR around 3C 249. 1 (iFu & Stocktonl2006l) . in 
the sense that both of them appear globally chaotic but locally 
ordered. 

3.2. Electron Density and Temperature 

Our IFUR datacube covers the [O II] AA3726,3729 doublet, 
and its spectral resolution (~1.8A) is high enough to resolve 
the two. We can determine the [O II] luminosity-weighted av- 
erage electro n densities (Ne\Q n]) from the intensity ratio of 
the doublet (Osterbrock 1989). Electron temperatures (T e ) 
can be measured using the [O III] (A4959+A5007)/A4363 ratio 
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FIG. 3. — [O III] A5007 radial velocity channel maps from the QSO-subtracted IFUR datacube. The central velocities, in km s~ x relative to that of the quasar 
NLR, are indicated in the lower-right corner of each panel. The circles show the extraction apertures of the various emission-line clouds discussed in the text. 
Some of the images have been allowed to wrap around to show both low surface brightness detail and high surface brightness peaks. The crosses indicate the 
position of the QSO before it was removed from the datacube. 



(Rom)- Most of the clouds in the IFUR field-of-view (FOV) 
have a low surface brightness, making it impossible to reliably 
measure the line ratios from a single 072 pixel. We therefore 
identified the peaks of the individual clouds, then binned pix- 
els within 0'.'3 of the peak. Figure [3] displays the extraction 
apertures for the eight resolved clouds. 

Figure |4j3 compares the extracted spectra from cloud e be- 
fore and after removing the QSO from the datacube (5 I2.U . It 
shows that the PSF subtraction technique can successfully re- 
move the QSO scattered light even for a cloud only Of! 6 from 
the quasar, as confirmed by the absence of the broad H/3 line 
in the residual spectrum. Note that since the data were taken 
under excellent seeing condition (~ Of! 6), PSF subtraction is 
only critical for the clouds very close to the quasar, namely 
a, e, and /. In addition, the line ratios for clouds a and / 
should be reliable even if the QSO removal were just barely 
correct, because the two are well separated from the quasar 
NLR in velocity. We are less confident about the line ratios 
of cloud e, since its velocity is almost centered on that of the 
NLR. However, the shape of the cloud as seen in the resid- 
ual datacube (see Fig.[3]as an example) suggests that the PSF 
removal cannot be too far off. Also at least the [O II] ratio 
should be trustworthy since the NLR emits very weakly in 
this doublet. 

To estimate the effect of using a fixed-size aperture on 
the covering fraction of the PSF (since the spatial resolu- 
tion gets better as the wavelength increases), we compared 
a QSO spectrum extracted from a d = 4'.'8 aperture with an- 
other one from a d = Of! 6 aperture 6 . A smoothed version of 
the ratio of the two was used to correct for the spectra of all 
the emission-line clouds. The atmospheric B-band absorption 
was corrected by dividing the spectra by an empirically deter- 
mined absorption law. The Galactic extinction on the line of 

6 Because the QSO is positioned at the upper left corner of the IFUR FOV, 
the extracted spectra were from the lower right quarter of the QSO only. We 
of course had to assume that the PSF was symmetric. 



sight to 4C 37.43 is A v = 0.072 (Schle gel et afl[l998h . Intrin- 
sic reddening due to dust associated with the cloud, whenever 
possible, was determined from the measured H7/H/3 ratio as- 
suming an intrinsic ratio of 0.46 8 (for case B recom bination at 
T = 10 4 K and N e = 10 2 cm' 3 : lOsterbrockll 1 9891) . However, 
there are four clouds, namely a, b, d, and g, where H7 is not 
securely detected in the IFU data. For clouds b and d, we used 
the reddening implied by the H7/H/3 ratio from the DEIMOS 
spectrum of SI (Fig. [I]). For a and g, we assumed an arbitrary 
reddening of Ay — 0.7. Both reddening eff ects were cor- 
rected using a standard Galactic reddening law (ICardelli et al.l 
fl989h . 

To measure the line ratios of the [O II] doublet, we used two 
Gaussian profiles constrained to be centered on the expected 
wavelengths of the doublet after applying the same redshift. 
As an example, Fig. @J? shows the profile of the [O II] dou- 
blet for the cloud e and the best two-Gaussian model fit. For 
the cases where there are two velocity components in a sin- 
gle spectrum, we first identify which component corresponds 
to the target cloud and which one is from contaminating light 
from other nearby clouds. Then we measure the velocity dif- 
ference between the two from the [O III] AA4959,5007 lines. 
Finally we freeze this parameter when fitting other lines, in- 
cluding the [Oil] doublet. For the DEIMOS spectra of Nl 
and SI, the [Oil] doublet is unresolved, so we measured the 
density-sensitive [S II] A6716/A6731 ratios. 

When there is a good [O III] A4363 line measurement, we 
derived T e and N e consistently using the IRAF routine TEM- 
DEN. Otherwise, a uniform T e = 10 4 K was assumed. 
N e [0 ll](A r e [S II]) would be approximately 15% (5%) higher, 
if T e = 15000 K. 

The physical properties of the EELR clouds are summa- 
rized in TableQ] The listed [O III] A5007 velocity dispersions 
have been corrected for the instrumental broadening (ctq ~ 35 
km s" 1 for IFUR; a w 67km s" 1 for DEIMOS). Reddening- 
free emission-line fluxes and 3-cr upperlimits are tabulated 
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FIG. 4. — Spectra of 4C 37.43 EELR e extracted from the IFUR datacube. (a) The He II A4686 to [O III] A5007 portion of the spectra. The upper and lower 
solid curves are the extracted spectra before and after removing the QSO from the datacube. The dashed curve shows the difference between the two, i.e., the 
scattered light from the QSO. For clarity, the dashed curve was shifted downward by 5 X 10 - erg cm" " 2 s" 1 A" 1 , (b) The [O II] AA3726,9 doublet of the 
EELR e. QSO scattered light has been removed. The dashed curves show the decomposition of the doublet, and the smooth solid curve is the sum of the two. 



in Table |2] as ratios to the H/3 flux. The quoted 1-a errors 
were derived from the covariance matrices associated with the 
Gaussian model fits, with the noise level estimated from line- 
free regions on either side of an emission line. For complete- 
ness, we have also included El in both tables. Most of the line 
fluxes were re-measured from the LRIS spect rum except for 
the [Q ui doublet ratio, which is quoted from lStockton et all 
d2002h . 

In general, the temperatures of the clouds are about 15000 
K, and the densities are a few hundreds cm~ 3 . The cloud e 
looks peculiar in this crowd: it has both the coolest tempera- 
ture and the lowest density, implying a pressure much lower 
than the average. The only other low-density cloud is Nl, 
where the [S II] ratio is at the low-density limit. 

3.3. Ionization Mechanisms 

Line ratio diagnostic diagrams are widely used to distin- 
guish different ionization mechanisms and to infer useful 
physical parameters of ionized clouds. Figure[5]presents three 
of these diagrams. Our line-ratio data are plotted against 
four ionization models, which were all computed by the 
multipurpose photoionization-shock code MAPPINGS: (1) 
a dusty radiation-pr essure dominated photoionization model 
dGroves et a l. 2004), (2) a two-phase photoionization model, 
(3) a s hock ionization model and (4) a "shock + precursor" 
model (iDopita & Sutherland! 19961) . 

These diagrams clearly show that both pure shock and 
"shock + precursor" are not as successful as the photoinoza- 
tion models in reproducing the observed spectra of the EELR. 
At first sight, Fig. [Sji might seem to imply that the "shock + 
precursor" model would possibly fit the data if we weight the 
precursor component more heavily. However, the high shock 
velocities (> 400 km s~ 4 ) required to reach the observed high 
[O lll]/H/3 ratios are inconsistent with the measured low ve- 
locity dispersions of those clouds (Fig.|2]i. In addition, mak- 
ing the precursor more luminous would only move the grids 
below where they are in Fig. [5J?, hence further off the data 
points. Therefore, we believe that it is unlikely that any of 
these clouds are dominated by shock ionization; instead, they 
are photoionized by the central QSO. In fact, it was argued 
before that at least El was photoionized by the QSO, as in- 
dicated by the low int ensity ratios of O IV A 1402/He II A1640 
and C iv A1549/He n dStockton et al.ll2002l) . 

Briefly, Model 1 simulates the photoionization of the sur- 
face layer (< a few pc) of a dense self-gravitating molecular 
cloud, where the hydrogen density near the ionization front is 



about 1000 cm -3 and the density structure is predominantly 
determined by radiation pressure forces exerted by the absorp- 
tion of photons by gas and dust. Detailed discussions of this 
model can be found in lGroves et al.l ( 120041) . 

The two-phas e photoionization m odel is similar to the 
Am ii model bv lBinette et alj dl996l) . It assumes that the ob- 
served emission-line spectrum is a combination of the emis- 
sion from two specie s of clouds. Guided by the results of 
IStockton etall d2002l) . we used the latest version of MAP- 
PINGS (3r) to model two isobaric components: (1) a high- 
excitation matter-bounded (MB) component with a density 
Nh ~ 1 cm~ 3 , a temperature T ~ 1.5xl0 4 K and an average 
ionization parameter U ~ 6 x 10~ 2 , and (2) a low-excitation 
ionization-bounded (IB) component with Nh ~ 400 cm -3 , 
T ~ 10 4 K and U ~ 2 x 10~ 4 . Unlike the A M/I model, 
shielding of the ionizing source by the MB component is not 
considered, i.e., both clouds are facing the same ionization 
field. Different combination ratios of the two components 
produce a sequence of spectra. We show in the diagrams a 
range of models having 20 to 80% of the total H(3 flux from 
the high-density IB component. 

Both photoionization models used simple power laws to 
represent the ionizing continuum (F v oc v a ). Model 1 as- 
sumed a single-index power- law from 5 to 1000 eV, and we 
show here two index values iGroves et al.l (f2()04) modelled, 
a = —1.2, —1.4. Model 2 uses a more realistic segmented 
power-law with a = —1.5 in the UV (8-185 eV) and a — 
- 1 .2 in the X-ray (185-5000 eV). The ionizing flux at the Ly- 
man limit is normaliz ed to 3.3 x 10 -17 erg cm -2 s _1 Hz -1 
dStockton et al.l l2002). These indexes are close to the canoni- 
cal value proposed for the typi cal spectral-energy distributio n 
of a QSO (i.e., a = -1.4; e.g- lFerland & Osterbrockl 19861) . 

We have run our two-phase model on two different sets of 
abundances: (a) th e same dust-deplet ed 1 Z & abundance set 
as used in Model 1 dGroves et alJ2004l); (b) the "standa rd" So- 
lar abundance set (Zq) of lAnders & Grevessd (Q989) scaled 
by -0.5 dex, following IStockton et all d2002l) . Note that the 
total gas metallicities in the two sets are both around one-third 
of Zq. Table [3] compares the two abundance sets. Sequences 
from both runs are shown in Fig. [5] In all diagrams, set b 
clearly provides a superior fit to the data than set a, meaning 
that the gas abundance in those clouds is more similar to Z Q . 
It also explains why there is a drastic offset 7 between the data 

7 This offset m ay imply that the Solar Ne/O abundance used by 
IGroves et aTU 2004) (log(Ne/0) = 0.61) is too high, by about a factor of three. 
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FIG. 5 . — Line ratio diagnostic diagrams distinguishing different ionization mechanisms: (a) [O II] AA3726, 9/H/3 vs. [O III] A5007/H/3, (b) He II A4686/H/3 
vs. [OlII] A4363/[OlII] A5007, and (c) [O II] AA3 726, 9/[Ne III] A38 69 vs. [O III] A5007/H/3. Overplotted are model grids from the dusty radiation-pressur e 
dominated photoionization model (black and grey; Groves et al. 2004), shock-only (blue) and "shock + precursor" models (green; IDopita & Sutherland] [T996T) . 
along with the two-component photoionization sequences with two different abundance sets (red lines with symbols). The dusty photoionization grids cover a 
range of ionization parameters (\og(U)) for metallicities of 1.0 and 2.0 Zq, and two power-law indexes for the ionizing continuum (black grids for a = —1.2, 
and grey ones for a = —1.4). In the shock models, a range of shock velocities (Vg [km s -1 ]) and magnetic parameters (B/n 1 / 2 [[i G cm -3 ' 2 ] ) are covered. 
The two-component photoionization sequences show a range of c omb ination ratios between the two components. The dashed and solid red lines represent the 
models adopting gas abundance set a and set b, respectively (see § l3.3l for details). The asterisks and diamonds, with increasing sizes, denote the positions where 
20, 50 and 80% of the total H/3 flux is from the high-density IB component. Line ratio measurements from the 4C 37.43 IFUR datacube are shown as blue open 
circles, and purple filled circles sh ow the measurements from either LRIS (El) or DEIMOS (Nl & SI) longslit spectra. For comparison, the four EELR clouds 
of 3C 249.1 (Fu & Stockton 2006) are plotted as red open squares. Arrows are reddening vectors. 

eled clouds. In the dusty model (Model 1) the hydrogen den- 
sity in a cloud varies from a few hundreds to ^1000 cm -3 . 
In this diagram, Model 1 cannot provide a good fit to the data 
unless an extremely low metallicity (< Q.2Z Q ) is used, which 
nonetheless contradicts the best fit grids in Fig.^ and Fig. [6] 
(§ 13.41) . However, the low-density MB component in the two- 
phase model can easily reproduce both the high temperature 
(~ 15000 K) and the low He II fluxes. Since most of the [O III] 
is emitted by the MB cloud, the influence from the IB cloud 
on temperature is negligible. We emphasize that the high- 
density IB component is indispensible to achieve the overall 
best fit, because the high [O II] fluxes cannot be explained by 
the high-excitation MB component alone. Note that since the 
emission lines involved in these diagrams are fairly insensitive 
to density effects, lowering density in the dusty model has the 
same effect as increasing the ionization parameter, and hence 
would not improve the fit to the data. 

The failure of Model 1 shows that these clouds are proba- 
bly not as de nse as a giant molecu lar cloud (GMC). In fact, 
as argued by Stockt on et alJ ((2002), the general lack of cor- 
relation between the morphology of the EELR and that of 
the stars in the host galaxy argues against them being GMCs, 
since GMCs are supposedly dense enough that their trajecto- 
ries should follow those of the stars and not be affected seri- 
ously by hydrodynamic interactions. 

The success of the simple two-phase model suggests that 
these clouds consist of a general low-density medium together 
with a clumpy distribution of much denser gas. This picture 
is also supported by the amount of reddening and the density- 
sensitive line ratios observed in these clouds (Table [TJ. An 
intrinsic reddening of Ay ~ 0.6 mag is generally seen in 
such clouds, which implies a hydrogen column density of 
~ 10 21 cm~ 2 . Assuming that the depth along the line-of- 
sight of the clouds is comparable to the resolved width of El 
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and the grids from Model 1 in Fig. |5j;. Therefore, hereafter 
the two-phase model (or Model 2) only refers to the run using 
abundance set b. 

Figure [5jz shows that both photoionization models can re- 
produce the strong oxygen lines, since such line ratios are 
mainly determined by the average ionization parameter. But 
the two begin to disagree in Fig.[2>, where weaker lines like 
He II A4686 and [O III] A4363 are involved. This discrepancy 
arises because of the different density structures of the mod- 



In fact, new Solar wind measurements i Gloeckler & Geiss 2007) have up- 
dated the Solar Ne/O with a new value (log(Ne/0) = 1.12), which is almost 
exactly three times lower. 
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FIG. 6.— Metallicity sensitive diagnostic diagrams — (a) [N II] A6584/Ha vs. [O III] A5007/H/3, (b) [N II] A6584/[0 II] AA3726, 9 vs. 
[O III] A5007/[O II] AA3726, 9. Keys are the same as in Fig. \5\ The photoionization girds here display a smaller range of ionization parameters but a wider 
range of metallicities (0.5, 1.0 and 2.0 Zq). In panel a the dashed black curve is the empirical starbursts/AGN dividing line i Kewley & Dopita 2002), and in 
panel b the dashed black lines are the empirical cuts used by Groves et al. (2006) to select the lowest metallicity Seyfert 2s (Z ~ IZq). In both panels, the top 
regions bordered by the dashed black curve/lines are where the low-metallicity Seyfert 2s are located. 



(1" ~ 5 kpc), the resulting number density is ^0.1 cm -3 . 
But we constantly measure a few hundreds cm -3 from either 
[O II] or [S II] ratios. A natural solution to this apparent con- 
tradiction is that clumpy dense clouds coexist with a much 
more diffuse medium. The former produces most of the low- 
ionization emission lines such as the [O II] AA3726,3729 and 
[S II] AA6717,6731 doublets, leading to the measured high 
de nsities; but the latter h olds most of the mass. As pointed out 
by lStockton et al.1 (120021) . the dense clouds need to be contin- 
uously supplied, possibly through compression of the diffuse 
medium by low-speed shocks, because (1) the two kinds of 
clouds are not in pressure equilibrium since the temperatures 
of the two are both around 10 4 K, and (2) the denser clouds 
are geometrically thin (~0. 1 pc, at least for the ionized part), 
thus cannot be stable on a timescale longer than 10 4 yr, given 
a typical sound speed of ~ 17 km s^ 1 . We will dis cuss the 
origin of the high-density medium in more detail in § 14.21 

3.4. Metallicity 

Diagrams involving the metallicity sensitive [N II] A6584 
line are shown in Figure [6] The only clouds where we have 
spectral coverage reaching the Ha-[N II] region are those ob- 
served with DEIMOS (Nl and SI). Although the two clouds 
are about 40 kpc apart in projected distance, their line ra- 
tios fall very close in Fig. [6] suggesting similar metallicities. 
iBoroson & Okd (1 19841) presented emission-line fluxes (their 
Table 2) of an El spectrum that covered the Ha-[N II] region. 
After dereddening, their data are consistent with the line ra- 
tios we measure from the LRIS spectrum of El. Because of 
their limited spectral resolution, they only listed the total flux 
from Ha + [Nil] AA6548,6584. We roughly decomposed the 
Ha and [N n] A6584 fluxes by assuming the case B value of 
HafH/3 = 2.9 and a theoretical ratio of the [N II] doublet of 
3.0. The line ratios of this region are also plotted in Fig.[6](la- 
beled as "Elp" to distinguish from the El region represented 
by the LRIS spectrum). 

The fact that different photoionization models with a same 
abundance 8 fall precisely onto one another in these plots 

8 Note that the total gas metallicities in abundance sets a and b (§ 13.31 used 



shows that such line ratios are almost entirely dependent on 
the abundance. Therefore, the grids from the dusty model 
(Model 1) can be used to measure the metallicities of the 
EELR clouds, even though it cannot reproduce some of the 
line ratios (refer to § 13.31 . 

By interpolating the grids of Model 1 in both panels of 
Fig. we estimated that 12 + log(0/H) = 8.55, 8.58, 8.42 
for Sl,Nl&Elp, respectively (Hereafter we will refer to 12 
+ log(0/H) as (O/H)). Since the metallicity determination also 
slightly depends on the density and Model 1 used a density of 
^1000, we have corrected the above values using the densities 
listed in Table Q] and the appr oximate (O/H) oc log(A^) scal - 
ing given in the Eqn. [4] of IStorchi-Bergmann et alJ (j 1998). 
The metallicities after the correction are (O/H) = 8.62, 8.68, 
8.47 for SI, Nl & Elp, respectively. The lower metallicity of 
Elp indicates that there might be some metallicity variations 
w ithin an EELR. 

IStorc hi-Berg mann et al.l d 19981) obtained two abundance 
calibrations for AGN NLRs from photoionization models. 
The calibrations involve similar line ratios as those in Fig. [6] 
and were shown to yield consistent abundance values with 
those from Hll regions for a sample of Seyferts. There- 
fore it is reassuring that the abundances obtained from these 
two calibrations are consistent with our simple interpolated 
results. The average o f the tw o abundance calibrations of 
IStorchi-Bergmann et alJ (119981) (their Eqs. [2], [3] and [4]) 
gives (O/H) ~ 8.64, 8.68 and 8.54 for SI, Nl and Elp, re- 
spectively. 

To compare the data with the Seyferts at similar redshifts 
from the Sloan Digital Sky Survey (SDSS), we delineated the 
regions where the lowest metallicity S eyfert 2s were fou nd 
in SDSS spectroscopic galaxy sample dGroves et alj|2006l) in 
Fig. [6] The regions are on the left side of the main Seyfert 
branch in both panels. It is clear that the three EELR clouds 
around 4C 37.43 all show significantly lower metallicities 
than typical NLRs of Seyferts ((O/H) > 8.8; iGroves et al] 

in the two-phase model are both about 1/3 Zq , and specifically, the nitrogen 
abundances are log(N/H) = —4.42 and —4.45, respectively. 
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4. DISCUSSION 
4. 1 . 77ie Low Metallicity of 4C 37.43 

Quasars are usually found in metal-rich environments 
((O/H) > (0/H)q = 8.93) (e.g., lHamann et al.1l2002h . The 
metallicity does not show a correlation with redshift, but 
seem s to increase with black hole (BH) masses (War neret al.1 
[2003b. The BH mass in 4C 37.43 is ab out 10 9 M w 
dLabita et al1l2006l) . which predicts (O/H) = 9.8 dWarner et al.1 
120031) . However, the low flux ratios of N v A1240/C iv A1549 
(~ 0.04) and Nv/HellA1640 (~ 0. 24) observed in the 
broad -line region (BLR) of this quasar (Kurasz kiewicz et all 
120021) both indicate a significantly lower metallicity, (O/H) 
« 8.4 {i.e., about 20 times lower than a typical quasar with 
Mbh = 10 9 M/t)) ac cording to the theoretical predictions of 
lHamann et al.1 {2002). Nevertheless, this BLR met allici ty is 
close to those measured in the surrounding EELR (§ I3.41 |, 

This low metallicity is unexpected for a quasar with a BH 
mass of 10 9 M if (1) the bulge-mass — BH-mass and the 
bulge-mass — metallicity correlations hold for this object and 
(2) the gas in the BLR and the EELR is from the ISM of the 
host galaxy of the quasar. There is no indication that the host 
galaxy of 4C 37.43 is particularly low mass: it is easily visi - 
ble in groun d-based and HST im aging ([Stockton et al. 2002), 
and, in fact, lLabita et"aT] d2006l) find that the BH mass for 
4C 37.43 from the host galaxy luminosity falls nicely between 
BH masses from the broad emission line widths calculated un- 
der two different assumptions regarding the geometry of the 
BLR (see their Table 4). 

Thus, whatever mechanism establishes the bulge-mass — 
BH-mass correlation likely was in operation in 4C 37.43 long 
before the current episode of quasar activity, and the BH al- 
most certainly underwent most of its growth phase early on. 
The low metallicity of the gas in the BLR and EELR therefore 
indicates that this gas has come from an external source, most 
likely a late-type gas-rich galaxy that has recently merged. If 
this is true, it provides the first direct observational evidence 
that gas from a merger can find its way down to the very cen- 
ter of the merger remnant, to within fueling distance of the 
active nucleus, although there has been a compelling statis- 
tical argument that this is the case for at least some QSOs 
dCanalizo & Stocktonll2001l) . 

4.2. The Origin of the High-Density Medium 

The viscosity of a fully ionized hydrogen plasma at N e = 1 
cm" 3 and T = 10 4 K is p w 9.5 x 10~ 7 g s" 1 cm" 1 dSpitzeri 
1962). Such a low viscosity implies a large Reynolds number 
if the scale and the speed of the flow are large. In fact, for 
a typical EELR cloud, R = Vm p Nl / p « 10 11 , given a flow 
velocity V = 200 km s -1 , a number density N = 1 cm -3 and 
a typical length scale I = 1 kpc. This means that turbulent 
motions can be easily triggered if the clouds are moving, con- 
verting part of its bulk kinetic energy into turbulent energy. 

The observed line-widths are around a ~ 50 km s _1 , while 
the sound speed in an ionized hydrogen plasma at 10 4 K is 
only c s ~ 17 km s — 1 , Shocks with a velocity ^50 km s _1 
must be common inside the cloud. If the pre-shock gas is 10 4 
K, an isothermal shock with a speed of 50 km s _1 is able to 
compress the gas by only a factor of 10 (pij p\ = M 2 , where 
M is the Mach number). But if the pre-shock gas is shielded 
from the ionizing source by the compressed post-shock gas, 
then the temperature of the pre-shock gas is probably ^100 



K, which corresponds to a sound speed of just ^2 km s . 
The same shock speed thus could increase the gas density by 
over 600 times, sufficient to explain the high density contrast 
(^300) between the two phases, as required by our best fit 
photoionization model (§ 13.3b . The total mass in the dense 
medium can stay relatively stable as long as the mass loss rate 
due to thermal expansion is comparable to the mass accretion 
rate through shock compression. We note that the amount of 
neutral gas due to this self-shielding should be inappreciable 
compared to the total mass of the cloud, since the filling factor 
of the dense materia l is known to be minuscule (<fi ~ 10~ 5 ; 
IStockton et al.l l2002). Also, for shocks at such a low speed, 
the amount of U V photons produced by the shock itself cannot 
significantly ionize the pre-shock neutra l gas (U ~ 6 x 10~ 5 
for Vs = 50 km s _1 and N = 1 cm~ 3 : lDopita & Sutherland! 
[1991 . 

Supersonic turbulence decays rapidly, typically over a 
crossing timescale. However, since these clouds are large, 
the crossing time is in fact on the same order-of-magnitude 
as the dynamical timescale of the EELR, r cross w 20(Z / 1 
kpc)/(u t /50 km s _1 ) Myr; so the picture of dense regions 
produced by shocks remains consistent with the timescales 
involved. 

4.3. Mass of the Ionized Gas 

The H/3 luminosity of an ionized cloud is proportional to 
N e N p V, where V is the volume occupied by the emitting ma- 
terial. So, if we assume Njj = N e = N p , then the total hydro- 
gen mass can be derived once we have an estimate of N e : M 
= A%m p f w dl/{ctKpN e hv) = 3.2 x 10 8 fBf>,-ie/N B)1 M , 
where m p is the proton mass, <Il the luminosity distance, an/3 
the effective recombination coefficient of H/3, hv the energy 
of a H/3 photon, /h/3,-16 the H/3 flux in units of 10~ 16 erg 
cm~ 2 s _1 and N Bi i the electron density in units of 1 cm -3 . 
We assu me a nominal density of 1 cm -3 because, as argued 
in § 13.31 most of the mass is in the low-density medium with a 
density of ~1 cm~ 3 even though the electron densities mea- 
sured from low-excitation lines are around hundreds cm -3 . 
The total mass of the entire EELR is ~ 3 x 10 10 M Q , given 
that the total H/3 flux is 1 x 10~ 14 erg cm~ 2 s -1 . 

One can also infer the cloud mass from the amount of dust 
observed and a given dust-to-gas ratio. The dust column den- 
sity is proportional to the amount of reddening we see along 
the line-of-sight. Assuming the Galactic dust-to-gas ratio, 
we obtain M = 5.8 x 10 21 cm~ 2 mag -1 EB-vm p d\9 2 = 
4.0 x 10 s Ay 9 2 M Q , where Ay is the intrinsic reddening 
in magnitude, 2 is the solid angle subtended by the cloud 
in □". This result is consistent with what we obtained from 
the H/3 luminosity if N e ~ 1 cm -3 , validating the existence 
of the low-density medium. Taking this one step further, if 
we require the two masses to be consistent, then it is easy to 
show that N e = 0.8 (/ H/ 3 _i 6 /0 2 )/Ay cm" 3 . The H/3 surface 
brightness of the EELR peaks at ~ 5 x 10~ 16 erg cm~ 2 s _1 
arcsec~ 2 for the brightest regions and at ~ 1 x 10 -16 erg 
cm~ 2 s _1 arcsec~ 2 for more average clouds (Fig. [2 given an 
[O lll]/H/3 ratio of 9). If we increase this value by a factor of 
2 to account for the smearing due to seeing, we end up with 
N e < 10 cm~ 3 for the brightest clouds and N e <2 cm~ 3 for 
the less luminous ones, given Ay = 0.8 mag. 

The mass of the EELR is on the same order of magnitude 
as that of the total ISM in the Galaxy, which implies that per- 
haps almost the entire ISM of the quasar host galaxy (or, more 
likely, that of the merging partner) is being ejected as well as 
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having been almost fully ionized by the quasar nucleus. 

4.4. The Driving Source of the Outflow 

The mass of the EELR, interpreted as an outflow, implies 
an enormous amount of energy. The bulk kinetic energy 
of the entire EELR is approximately Eke — MV 2 /2 = 
6.2 x 10 5 M10V250 ergs, where M\o is the total mass in 
units of 10 10 Mq and V 250 = V/250 km s _1 . We used 250 
km s" 1 as the nominal velocity scale because most of the 
mass is in the brightest two concentrations, which are both 
blueshifted by about 250 km s" 1 . The kinetic energy of the 
unresolved kinematic substructures ("turbulent" energy) can 
be estimated from the measured line widths, E tur — Ma 2 = 
5.0 x 10 56 Mi cr| ergs (ct 50 = a/50 km s" 1 ). The emission- 
line gas also has a thermal energy of Eth = 1.5MkT/m p = 
2.5 x 10 55 Mi T 4 ergs (T 4 = T(K)/10 i K), which is just 
about 0.4% of the kinetic energy. Similarly, the total momen- 
tum is p = MV = 5 x 10 50 MioV 2 50 dyne s. Whatever the 
driving source of this outflow is, it must have deposited this 
enormous amount of energy and momentum into the EELR in 
a short time (a few x 10 7 yr, the dynamical timescale of the 
EELR). Such high i nput rates are consis tent with an outflow 
driven by a quasar dFu & Stocktodl2006l) . 

Although in most cases the EELR morphology seems inde- 
pendent of the structure of the radio source, the fact that both 
the occurrence and l uminosity of the EELR increa se with the 
radio spectral index (Sto ckton & MacKentvl fl987) suggests a 
link between the extended emission and the radio jet. Re- 
cent simulations of jet-cloud interactions in elliptical galaxies 
show that wide-solid-angle "bubbles" could accompany the 
launching of the jet (M. Dopita 2007, private communication). 
This is an intriguing picture, because such bubbles provide a 
natural way to expel clouds from the vicinity of the quasar to 
radii of a few tens of kpc without resulting in any morpho- 
logical similarities between the expelled material and the jet, 
while also producing the observed connection between radio 
outflows and EELRs. 

This mechanism for globally expelling a quantity of gas 
comparable to the total ISM of a reasonably massive galaxy 
is certainly of interest in the context of current speculation re- 
garding quasar winds as a means of initially establishing the 
observed bulge-mass — black-hole-mass correlation (regard- 
less of the extent to which it has to be maintained by heat- 
ing the surrounding medium via less violent nuclear activity). 
The main differences are that, (1) in 4C 37.43, the black-hole 
has already achieved a mass of ~ 10 9 M Q , presumably in a 
previous quasar phase, and it is currently ejecting rather low- 
metallicity gas that likely comes from a late-type, gas-rich 
merging companion; and (2) the ejection of the gas is a direct 
consequence of large-scale shocks produced by the radio jet, 
rather than due to radiative coupling of the quasar luminosity 
to the gas. 

5. SUMMARY 

Our results suggest the following overall picture for the ori- 
gin of the EELR around 4C 37.43. A large, late-type galaxy 
with a mass of a few x 10 10 M Q of low-metallicity gas has 
recently merged with the gas-poor host galaxy of 4C 37.43, 
which has a ~ 10 9 M Q BH at its center. The low metallicity 
of the BLR gas indicates that gas from the merging compan- 
ion has been driven to the center during the merger, trigger- 
ing the quasar activity, including the production of FR II radio 
jets. The initiation of the jets also produces a wide-solid-angle 



blast wave that sweeps most of the gas from the encounter out 
of the galaxy. This gas is photoionized by UV radiation from 
the quasar, and turbulent shocks produce high-density (~ 400 
cm -3 ) filaments or sheets in the otherwise low-density (~ 1 
cm -3 ) ionized medium. The EELR will have a lifetime on 
the order of 10 Myr, i.e., comparable to that of the extended 
radio source. 

In more detail, we can summarize the major conclusions of 
this paper as follows: 

1. The EELR of 4C 37.43 exhibits rather complex kine- 
matics which cannot be explained globally by a simple 
dynamical model. 

2. The [O II] or [S II] electron densities of the clouds range 
from 600 cm" 3 to less than 100 cm -3 . The Rom 
temperatures are mostly ~ 1.5 x 10 4 K. The cloud 
(e) having the lowest temperature (T ~ 1.2 x 10 4 K) 
also shows the lowest density, indicating a lower-than- 
average pressure. 

3. The spectra from the clouds are inconsistent with shock 
or "shock + precursor" ionization models, but they are 
consistent with photoionization by the quasar nucleus. 

4. The best-fit photoionization model requires a two-phase 
medium, consisting of a matter-bounded diffuse com- 
ponent with a unity filling-factor (N ~ 1 cm" 3 , T ~ 
15000 K), in which are embedded small, dense clouds 
(N ~ 400 cm" 3 , T ~ 10 4 K), which are likely con- 
stantly being regenerated through compression of the 
diffuse medium by low-speed shocks. 

5. The metallicity of the EELR ((O/H) < 8.7 ) is similar to 
that o f the rare low-metallicity Seyferts dGroves et al.l 
2006), as implied by [N II] A6584 line ratios and the 
overall best-fit photoionization model. Previous results 
show that the BLR has a similarly low metallicity, in- 
dicating a common (external) origin for the gas, which 
also presumably fuels the current quasar activity. 

6. The photoionization model gives a total mass for the 
ionized gas of 3 x 10 10 M Q . The total kinetic energy 
implied by this mass and the observed velocity field is 
- 2 x 10 58 ergs. 

7. The strong correlation of luminous EELRs with steep- 
spectrum radio-loud quasars, coupled with the general 
lack of significant correlation between the EELR and 
radio morphologies implies that the highly collimated 
radio jets are accompanied by a more nearly spherical 
blast wave. 

8. Since the mass of the ionized, apparently ejected gas in 
4C 37.43 is comparable to that of the ISM of a mod- 
erately large spiral, this object provides a local ana- 
log to (though likely not an example of) the hypoth- 
esized "quasar-mode" ejection that may be instrumen- 
tal in initially establishing the bulge-mass — black-hole- 
mass correlation at high redshifts. 



We thank Tracy Beck and Gelys Trancho for helping with 
the GMOS observations, Marc Kassis and Chuck Sorenson 
for supporting the DEIMOS observations, and Mike Dopita, 
Lisa Kewley, and Brent Groves for helpful discussions on the 



Extended Emission-Line Region of 4C 37.43 



11 



photoionization models. We thank Steve Rodney and Ryan 
Foley for kindly providing the DEIMOS standard star obser- 
vation files. We also thank the anonymous referee for a care- 
ful reading of the manuscript, for his or her concern that the 



Allington-Smith, J., et al. 2002, PASP, 1 14, 892 

Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 

Binette, L., Wilson, A. S., & Storchi-Bergmann, T. 1996, A&A, 312, 365 

Boroson, T. A., & Oke, J. B. 1984, ApJ, 281, 535 

Boroson, T. A., Persson, S. E., & Oke, J. B. 1985, ApJ, 293, 120 

Canalizo, G., & Stockton, A. 2001, ApJ, 555, 719 

Cappellari, M., & Copin, Y. 2003, MNRAS, 342, 345 

Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 

Crawford, C. S., Fabian, A. C, & Johnstone, R. M. 1988, MNRAS, 235, 183 

Crawford, C. S., & Vanderriest, C. 2000, MNRAS, 315, 433 

Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 604 

Dopita, M. A., & Sutherland, R. S. 1996, ApJS, 102, 161 

Durret, E, Pecontal, E., Petitjean, P., & Bergeron, J. 1994, A&A, 291, 392 

Faber, S. M., et al. 2003, Proc. SPIE, 4841, 1657 

Fabian, A. C, Crawford, C. S., Johnstone, R. M., & Thomas, P. A. 1987, 

MNRAS, 228, 963 
Ferland, G. J., & Osterbrock, D. E. 1986, ApJ, 300, 658 
Friaca, A. C. S., & Terlevich, R. J. 1998, MNRAS, 298, 399 
Fu, H., & Stockton, A. 2006, ApJ, 650, 80 
Gloeckler, G, Geiss, J. 2007, Space Science Reviews, in press. 
Groves, B. A., Dopita, M. A., & Sutherland, R. S. 2004, ApJS, 153, 9 
Groves, B. A., Heckman, T. M., & Kauffmann, G. 2006, MNRAS, 371, 1559 
Hamann, F, Korista, K. T., Ferland, G. J., Warner, C, & Baldwin, J. 2002, 

ApJ, 564, 592 

Hook, R., Lucy, L., Stockton, A., & Ridgway, S. 1994, Space Telescope 
European Coordinating Facility Newsletter, Volume 21, p. 16, 21, 16 



paper be presented in a way that would make it accessible to 
a wide audience, and for cogent comments to help us clarify 
and improve the presentation. This research has been partially 
supported by NSF grant AST03-07335. 



Hook, I. M., J0rgensen, I., Allington-Smith, J. R., Davies, R. L., Metcalfe, 
N., Murowinski, R. G., & Crampton, D. 2004, PASP, 116, 425 

Hopkins, P. F, Hernquist, L., Cox, T. J., Di Matteo, T., Robertson, B., & 
Springel, V. 2006, ApJS, 163, 1 

Kewley, L. J., & Dopita, M. A. 2002, ApJS, 142, 35 

Kuraszkiewicz, J. K., Green, P. J., Forster, K., Aldcroft, T. L., Evans, I. N., & 

Koratkar, A. 2002, ApJS, 143, 257 
Labita, M., Treves, A., Falomo, R., & Uslenghi, M. 2006, MNRAS, 373, 551 
Miller, P., Rawlings, S., & Saunders, R. 1993, MNRAS, 263, 425 
Oke, J. B., et al. 1995, PASP, 107, 375 

Osterbrock, D. E. 1989, Astrophysics of Gaseous Nebulae and Active 
Galactic Nuclei, Mill Valley, California: University Science Books 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 

Spitzer, L. 1962, Physics of Fully Ionized Gases, New York: Interscience (2nd 
edition), 1962 

Stockton, A. 1976, ApJ, 205, LI 13 

Stockton, A., Fu, H, Henry, J. P., & Canalizo, G., 2006a, ApJ, 638, 635 
Stockton, A., Fu, H, & Canalizo, G. 2006b, New Astronomy Review, 50, 694 
Stockton, A., & MacKenty, J. W. 1987, ApJ, 316, 584 

Stockton, A., MacKenty, J. W., Hu, E. M., & Kim, T.-S. 2002, ApJ, 572, 735 
Storchi-Bergmann, T, Schmitt, H. R., Calzetti, D., & Kinney, A. L. 1998, AJ, 
115,909 

Warner, C, Hamann, R, & Dietrich, M. 2003, ApJ, 596, 72 



TABLE 1 

Properties of 4C 37.43 EELR clouds 



R poinn 


v 






H/3 x 10 17 


in nl 

IKJ 11J 


rs nl 


in ml 

111J 


N b 

1 ' e 


T 

e 




(km s _1 ) 


(km s -1 ) 




(erg cm 2 s ^) 


3726/3729 


6717/6731 


(4959+5007)/4363 


( cm ~~ ^ ^ 


(K) 


a 


-363 


68 


(0.700) 


4.4 ±0.4 


1.17 ±0.22 






580 ± 300 




h 


-121 


133 


(0.542) 


6.9 ±0.3 


1.08 ±0.10 






460 ± 130 




c 


-150 


79 


0.540 


8.2 ±0.4 


0.86 ± 0.09 




77 ± 24 


230 ± 120 




d 


-55 


116 


(0.542) 


5.9 ±0.4 


1.07 ±0.24 






450 ± 310 


e 


-13 


69 


1.534 


164.6 ±2.2 


0.68 ±0.05 




107 ± 15 


20 ±50 


12500+™ 


f 


85 


79 


1.362 


23.1 ±0.7 


0.98 ±0.11 




69 ±26 


390 ± 160 


isooolf™ 


g 


167 


101 


(0.700) 


9.4 ±0.5 


0.89 ±0.12 






230 ± 140 


h 


214 


63 


1.380 


23.3 ±0.9 


1.04 ±0.11 




57 ± 14 


500 ± 170 


16400+?°°° 


El 


-250 


80 


0.655 


317.2 ±3.2 


0.97 ± 0.03 




75 ±4 


375 ± 40 


14500 ± 400 


SI 


-73 


162 


0.542 


6.8 ±0.4 




1.24 ±0.10 




190 ± 120 




Nl 


-59 


148 


0.237 


11.7±0.3 




1.46 ±0.11 


66 ±25 


< 60 


1 5300+ 4300 



a Intrinsic reddening. The values in parentheses are not directly measured from the H7/H/3 ratio, since the H7 lines in those spectra are not well- 
detected. See § 13.21 for details. b If there is a good measurement of the temperature sensitive [OlII] intensity ratio, the electron density is derived 
consistently with the Rom temperature. However, the quoted i-cr error bars do not include the uncertainties of the other parameter. In other cases, N e 
is derived assuming T e = 10 4 K. 
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TABLE 2 

Line Ratios of 4C 37.43 EELR Clouds Relative to H/3 



Region [NeV]A3426 [OII]A3726 [OII]A3729 [NeIII]A3869 [OIII]A4363 [HeII]A4686 [OIII]A5007 



Ha 



[NII]A6584 [SII]A6716 [SII]A6731 



a 




< 


1.51 


1 


53 ± 0.18 


1.31 ±0.19 





85 ±0 


27 




< 0.41 




< 0.29 


9 


30 ±0 


13 


b 




< 


0.98 


2 


35 ± 0.15 


2.17 ±0.14 


1 


04 ±0 


14 




< 0.28 




< 0.18 


6 


88 ±0 


07 


c 




< 


0.61 


2 


10 ±0.09 


2.43 ±0.23 


1 


27 ±0 


11 





15 ±0.05 




< 0.13 


8 


50 ±0 


13 


d 




< 


1.03 


1 


88 ± 0.28 


1.76 ±0.29 





97 ±0 


13 




< 0.30 




< 0.20 


8 


21 ±0 


08 


e 


1 


85 


±0.08 


1 


07 ±0.07 


1.57 ±0.07 


1 


14 ±0 


05 





15 ±0.02 





33 ±0.01 


12 


01 ±0 


01 


f 




< 


0.69 


2 


04 ±0.21 


2.09 ±0.08 





98 ±0 


09 





14 ±0.05 





18 ±0.03 


7 


37 ±0 


11 


g 




< 


0.62 


1 


67 ±0.20 


1.89 ±0.11 





80 ±0 


12 




< 0.23 





21 ±0.04 


9 


73 ±0 


10 


h 


1 


00 


±0.22 


1 


59 ± 0.12 


1.53 ±0.12 


1 


46 ±0 


11 





19 ±0.05 





19 ±0.03 


8 


30 ±0 


03 


El 





36 


±0.02 


1 


42 ± 0.03 


1.46 ±0.03 





96 ±0 


02 





19 ±0.02 





20 ±0.02 


10 


50 ±0 


06 


SI 




< 


0.38 




4.42 ± 0.10 a 





73 ±0 


07 




< 0.17 





22 ±0.05 


8 


74 ±0 


07 


Nl 




< 


0.14 




3.53 ± 0.05 a 





81 ±0 


01 





16 ± 0.06 





15 ± 0.03 


8 


01 ±0 


01 



o 
o 

o 

3 



3.16 ±0.03 
3.16 ±0.03 



1.50 ±0.03 
1.56 ±0.03 



0.55 ±0.03 0.44 ±0.03 
0.64 ±0.03 0.44 ±0.03 



a These values are lolal fluxes in the [O ill doublet, since the doublets are unresolved in these spectra. 
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TABLE 3 
Abundances of key elements 



Set 1 


H 


He 


N 


O 


Ne 


a 


0.00 


-0.99 


-4.42 


-3.53 


-3.92 


b 


0.00 


-1.01 


-4.45 


-3.57 


-4.41 



All abundances are logarithmic with respect to hydrogen. 



